Are Domain Walls in 2D Spin Glasses described by Stochastic Loewner Evolutions? 
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Domain walls for spin glasses are believed to be scale invariant; a stronger symmetry, conformal 
invariance, has the potential to hold. The statistics of zero-temperature Ising spin glass domain 
walls in two dimensions are used to test the hypothesis that these domain walls are described by a 
Schramm-Loewner evolution SLE K . Multiple tests are consistent with SLE K , where K = 2.32 ±0.08. 
Both conformal invariance and the domain Markov property are tested. The latter does not hold in 
small systems, but detailed numerical evidence suggests that it holds in the continuum limit. 



The geometrical characterization of physical objects is 
central to much of our understanding of their energetics 
and dynamics. The relevant geometries can be as simple 
as points or gently curved surfaces. Many objects are not 
well-described by an integer dimension, but have a scale- 
dependent measure that can be represented by a fractal 
dimension. For example, continuous phase transitions in 
homogeneous systems have nonanalytic behavior consis- 
tent with fractal dimensions for the surfaces that separate 
phases. Evidence for fractal domain walls is seen in scat- 
tering experiments and numerical simulations. In mod- 
els of glassy systems with quenched disorder (frozen-in 
random fields), analytic work and numerical simulations 
indicate that domain walls can be scale-invariant and 
fractal at low temperatures. In two-dimensional homo- 
geneous systems, the additional symmetry of conformal 
invariance often applies and yields detailed predictions 
for critical exponents, the effects of boundary conditions, 
and a background for physical explanations. 

The conjunction of conformal invariance with the pres- 
ence of a domain Markov property (DMP) in statistical 
mechanics models has led to an even more complete - 
and in several cases mathematically rigorous - descrip- 
tion of fractal curves such as loop-erased self-avoiding 
walks, percolation hulls, and domain walls at phase tran- 
sitions in two dimensions in the scaling (i.e., continuum) 
limit [J, 0]. Schramm showed that when both proper- 
ties are present the probability measure on these curves 
is described by a Schramm-Loewner evolution SLE K [3|. 
Random sequences of simple conformal maps can be used 
to generate the fractal curves with the correct measure, if 
the real-valued driving function that underlies the maps 
is a Brownian motion. The diffusion coefficient k of 
the Brownian motion uniquely parameterizes the process 
and is related to the fractal dimension of the curve via 
df = l+«/8. This deep connection has led to very precise 
characterization of these curves for pure systems such as 
g-states Potts model, 0(n) models and percolation. 

An outstanding question is whether SLE can be ap- 
plied to other systems. Numerical evidence has been pre- 
sented, for example, that SLE describes certain isolines 
in 2D turbulence The broad question of whether and 
how conformal invariance, a necessary condition for SLE, 



applies to disordered systems is still very much open. At- 
tempts to extend the apparatus of conformal field theory 
to systems with quenched disorder, a notably difficult 
subject [g|j have suggested some numerical tests, such 
as finite size scaling [6j]. A positive result was obtained 
recently for surface wavefunction multifractality at the 
2D localization transition with spin-orbit symmetry 0]. 
Most importantly for our work, Amoruso, Hartmann, 
Hastings and Moore [j| have suggested that domain walls 
in the 2D spin glass have conformal invariance. 

Here we directly investigate whether the domain wall 
statistics converges to an SLE K . We apply several tests. 
We examine the winding of domain wall around a cylin- 
der, as well as the angular distribution of the curves and 
the dipolar SLE hitting probability. We use an iterated 
slit map (discretized inverse Loewner evolution) to deter- 
mine the driving function and test whether it converges 
to Brownian motion. We directly test the DMP by com- 
paring precise domain wall statistics in "whole" and "cut" 
domains. To determine the significance of these tests, we 
carry out the same analyses for the loop-erased random 
walk (LERW), which has SLE2 as a scaling limit, and for 
paths on minimal spanning trees (MST), which are not 
conformally invariant [lfj. We find that, for one choice 
of boundary conditions (BCs), the spin glass domain wall 
passes all tests with a consistent value of n. 

We study the domain walls in a 2D Ising spin glass 
with Gaussian disorder. We use the Edwards- Anderson 
Hamiltonian Ti. = — Jij s i s ji where Sj = ±1 and 

the Jij are each chosen from a Gaussian distribution with 
zero mean. The glass transition is at T — 0; we study the 
minimum energystates at T = using an exact optimiza- 
tion algorithm [12| and sample over disorder realizations. 
There are two ground states, connected by a global spin 
flip, in any finite sample. In the scaling picture based on 
domain walls and droplets, introduced by McMillan 
Bray and Moore and Fisher and Huse 15 1, there are 
two ground states in the thermodynamic limit 0,0. 
Domain walls (DWs) separate these two ground states. 
The domain wall energy scales as Edw ~ L e flU 
for DWs defined at scale L, with 6 = -0.28(1). 

We work on a triangular lattice that has W spins in 
each of L rows, as sketched in Fig. [TJ Our samples 
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are cylindrical, with periodic rows. One can uniquely 
describe ground state pairs by the bond satisfactions 
aij = sga(JijSiSj), where n A c7y = sgn(n A J y ) for any 
elementary triangle A. Periodic BCs result from fixing 
II(y') e p a o"ij = Huj) g Po ,sgn( Jij), where p a are the sets of 
boundary bonds on the upper (a — 1) and lower (a = 2) 
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gives antiperiodic BCs (equivalent to a change in the 
sign of the Jij along a column of bonds). Comparing 
ground states for periodic and antiperiodic BCs gives 
a domain wall: a simple path on the dual lattice that 
crosses bonds whose satisfaction differs between the two 
BCs. A domain wall is the minimizer of the cost func- 
tion 2 J2(ij)e-y Fi J '' with F v = Jij s i s j ar, d s° the spins 
in the periodic ground state, over open paths 7 from the 
bottom to the top of the cylinder. In the ground state, 
the cost of any closed loop is positive fl8| . 

We refer to the domain wall found using this par- 
ticular periodic-antiperiodic BC comparison as "float- 
ing" (F-PA), as the endpoints of the domain wall are 
not fixed. We also consider a periodic-antiperiodic BC 
change where the domain wall at one end is locally con- 
strained to a single chosen bond on the lower boundary 
(L-PA), i.e., a given o; L j is reversed on the lower boundary. 

The choice of L and W give the cylinder shape, with 
the circumference given by X — W and the length by 
Y = (V3/2)L. We find that the averages converge for 
W > 4L. The results for the first quarter of the path for 
W = L agree with those for W ^ L, within our accu- 
racy. For comparison, we also study LERW curves and 
paths between two points in the MST, both on honey- 
comb lattices; LERW curves have dimension df = 5/4 
in the continuum limit and MST paths appear to have a 
fractal dimension df = 1.217(3) [9(|. 

We estimate df for the domain wall by computing the 
mean total path length S(L) of the domain wall, com- 
paring with S(L) ~ L d f , the overline indicating aver- 
ages over w 10 4 samples at sizes up 5 x 10 5 spins, and 
also by computing the sample averaged distance from 
the origin as a function of partial path length. We find 
dj_= 1.28(1), in statistical agreement with previous work 
[lH, for both F-PA and L-PA BCs. 

We test conformal invariance and consistency with the 
SLE description by measuring the winding of the F-PA 
domain wall around a long cylinders with Y ^> X. The 
prediction from SLE is that the variance of the transverse 
displacement x of the end point from the starting location 
is (x 2 ) = £ (df - l)XY. We studied cylinders with 8 < 
W < 32 and up to L = 800 for at least 10 4 samples at 
each of at least eight values of L. Our data is consistent 
with \ {x 2 ) linear in XY, with a coefficient of 0.27± 0.01 
for L > AW, in agreement with conformal invariance, 
again for both F-PA and L-PA BCs. 

A powerful result from SLE is a prediction of the 
probability that a curve generated by SLE will pass to the 
left of a given point at polar coordinates (R, 4>) (see Fig. 
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L=W=400, R=50, F-AP 
L=W=720, S=90, F-AP 
L=W=720, R-3m, F-AP 
L=128, W=1024,fl=64, 
L-AP (fixed end) 
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Figure 1: Plots of P(4>, R)— P2 ((/>), where P(<j>, R) is the proba- 
bility for the domain wall to pass to the left of a point with po- 
lar coordinates (R, <f>) (see inset). The magnitude of statistical 
errors (not shown) is consistent with the apparent fluctuations 
of the data lines. The predicted P2(<t>) = [4>— | sin(2(/j)]/vT for 
k = 2 is subtracted to display small variations clearly. The 
data from F-PA paths agrees with SLE predictions for k in the 
range 2.24 < k < 2.40, while L-PA paths give k » 2.85(10). 
Inset: A domain wall of length S = 9 in a sample with L — 4 
rows and W = 6 columns. 



Q]for notation). Given scale invariance, the probability 
that the curve passes to the left of (i?, </>) depends only 
on (j>, and the theory of SLE can be used to predict (l9j 
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where F12 is the hypergeometric function and K is the 
diffusion parameter from continuum SLE. Our results 
(Fig. [1]) for P((p) depend on the choice of BC. For F- 
PA BCs the measured P(4>) is most consistent with the 
analytical form in the range k^ s = 2.32 ±0.08, consistent 
with the relation n = 8(df — 1). For L-PA domain walls 
measured from the fixed end, we find n^g = 2.85 ± 0.10. 
We find the same two BC-dependent values [13] using an- 
other test: a comparison of the distribution for the dis- 
placement between the DW endpoints with the form pre- 
dicted using dipolar SLE 



2 11 ] . which describes the limit 



X/Y —> 00 for SLE curves that start at a given point a 
and terminate on the upper boundary. Constraining the 
domain wall to start at a given point (L-PA BCs) rather 
than choosing domain walls that start at a point (i.e., 
conditioning on a) with F-PA BCs changes the effective 
k. Domain walls with fixed endpoints are not consistently 
described by SLE K over their entire length. 

The boundary conditions appear not to affect the frac- 
tal dimension, but clearly do affect more subtle aspects 
of the geometry. A similar result holds for the LERW: 
absorbing and reflecting boundary conditions both are 
consistent with df = § [2|, but P(4>) is well fit by Eq. (U) 
with k = 2 for absorbing BCs, as expected, in contrast 
with P(4>) ~ 4>/ir for reflecting BCs. 
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Note that the L-PA DW energy approaches a constant 
as L — ¥ oo, in contrast with £bw ~ L 6 for F-PA DWs. 
This difference holds in general for 8 < 0, as can be seen 
by summing L-PA domain wall energies defined over ge- 
ometrically increasing scales connecting the localized re- 
gion to the large scale. Essentially, the L-PA constraint 
gives an O(l) correction to £bw from the cost of the 
single bond at the localized end. Apparently, the opti- 
mization over successive scales distorts the curve from 
the form expected from SLE, while optimization over a 
single global scale gives results consistent with SLE. 

To more carefully inspect the correspondence with 
SLE, we have used a discrete Loewner evolution to map 
the domain walls, represented by sequences of points 
zf = x" +iyi, i = 1 . . . S, in the complex half-plane, onto 
a real- valued sequence defined at discrete U and 

studied the sample statistics and correlations of the in- 
terpolated £(t). For continuous curves generated by SLE, 
the underlying function £ (t) is Brownian motion with dif- 
fusion constant K. The sequence is initialized by setting 
t = and f(io) = 0. We then recursively map the 



sequence {z\ 



} to the transformed and short- 



ened sequence {z l i+1 , . . . , z l s } using the map appropriate 
for dipolar SLE (first defining A, = %y\/2Y), 



U = h-i -2(y/ 7 r) 2 log[cos(A I )] ; £(*i) = 



(2) 



2L 



Zj = H cosh 1 <J cosh 



?r(z-&) 

2L 



I cos(Ai; 



These maps are a sequence of slit maps that successively 
remove the first point from the sequence (see Fig. [21) and 
maintain the hydrodynamic normalization used in SLE. 

The simplest test for the diffusive property of is 
to examine its distribution at fixed times. Our data for 
L 2 /5 > t > 50 are consistent with a Gaussian distribu- 
tion for with variance £ 2 (£) = n e fft (Fig. [2]). We 
have confirmed that higher cumulants satisfy £ 2n (i) = 

(2n)\\ (j; 2 {t)J for n = 2,3,4, within numerical error, 

for the same range of t. For computations of £(t) that 
start from a free end of a domain wall (F-PA boundary 
conditions or the free end of L-PA BCs), we find nf s = 
2.24 ± 0.08, while for computations starting from the lo- 
calized end with L-PA BCs, we estimate k^ s « 2.85±0.1, 
consistent with our estimates from P(<fr, R). We note that 
£ 2 (i) is also nearly linear in time for paths on the MST, 
even though such paths are not conformally invariant, 
but the coefficient is not consistent with the fractal di- 
mension (see [13] for MST winding angle results). 

We have also tested the Markovian property for 
i.e., that the changes in £(t) depend only on the cur- 
rent value of and not on previous values. We 
studied the correlation function Cd(n) = — 
£(ii+n)][£(^+l) — £(tt)]) at intermediate times; it decays 
rapidly (by a factor of » 100 over the range n = 2 to 
n = 8) for both the spin glass and for the LERW [IJ. 
Note that there must be short term correlations in £(t) 
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L=256, from constrained end 
L=128, from constrained end 
t=128, floating BCs 
L=64, floating BCs 



Figure 2: P lot of an effective diffusion constant K e « = 
£ 2 (2t)-£ 2 (t)/t, for W > 4L. Lines indicate K = 2.24, 2.32, 
2.40, 2.85, and 3.00. The range 2.24 < « < 2.40 fits the data 
for curves with F-AP BCs, while 2.85 < K < 3.00 describes 
the diffusion measured from a constrained domain wall end. 
Inset: Part of a sample conversion of a domain wall in the 
2D Ising spin glass to a sequence £(fi), i ~ ^ ■ ■ ■ S. The left 
curve is the initial domain wall with £(0) = 0, while the red 
[lighter] curve is the remainder after 500 applications of the 
dipolar map, giving £(tsoo ~ 7239.4) » 101.5. 



on the lattice, as there are forbidden sequences of "turns" 
for the domain wall. 

Given that the Ising spin glass DW passes several SLE 
tests one must examine the domain Markov property 
(DMP) in a disordered system. Let us call Poilab) the 
probability that the DW happens to coincide with the 
curve 7ab in a domain D (where a,b are two given bound- 
ary points). The DMP [22] states that if one conditions 
this probability on a piece j a c of the curve, then the 
probability for the rest of the curve j c b is identical to the 
original probability on the cut domain D\j ac conditioned 
on curves starting at c, i.e: 



PD{lcb\lac) = PD\ lac hcb\c) 



(3) 



Cutting the domain removes bonds that cross the seg- 
ment 7 ac . In pure statistical systems Eq.[3]is an identity, 
given proper BCs. One can easily check that the DMP 
holds in a single realization of disorder. However this 
property does not survive disorder averaging (as condi- 
tioned probabilities are ratio of probabilities) except for 
percolation SLE K= 6 (because of locality). 

To evaluate the deviations from the DMP we have com- 
puted numerically the ratio of sums of the two proba- 
bilities in Eq. (J3j) - We generate domain walls in both 
whole cylinders and in cylinders cut by all paths 71 of 
a chosen path length Si, We sample at least 3 x 10 7 
disorder configurations to estimate the ratio r(7i,7a) = 
S 7bcD72 P D\ 7l (7fcc|c)/E 7bcD72 ^u(76c|7i). with 71 start- 
ing at the lower boundary, 72 a subpath of 7h c of path 
length S2, and 71 connecting to 72 at c. If DMP holds 
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Figure 3: Plot of C(x), cumulative probability of ranked val- 
ues for r(7i, 72), as defined in the text. Large deviations from 
r = 1, as clearly seen for R-LERW, indicate a failure of the 
domain Markov property. 



strictly, ^(71, 72) = 1. We summarize our data in Fig. 
O where we plot the cumulative probability C{x) = 
E 7 i, 72 |r( 7l ,7 2 ) <a; p d(ii, 72) that r is less than x. We find 
that in small samples with L — 8, r is statistically dis- 
tinct from unity for larger I72I. The largest deviations are 
seen for 72 near to and parallel to 71. For comparison, 
we show results of the same analysis for the LERW with 
both absorbing BCs (A-LERW), where r is unity within 
statistical error, and for reflecting BCs (R-LERW), where 
r clearly deviates from unity. In ISG simulations, C{x) 
is quite close to the curve for A-LERW. Our data cannot 
rule out the possibility of DMP holding in the continuum 
limit. 

It is tempting to conjecture that the emergence of the 
DMP in the continuum limit follows from the existence 
of "principal" minimizers, separated from each other on 
a scale L. These are the basins of attraction for min- 
imizing paths: if the start of a DW is displaced from 
the minimizer's start by a scale t < L, the DW merges 
with principal minimizer within a distance of order I [1 61 ] - 
These minimizers are a result of finding shortest paths in 
a graph with negative weights (but no negative weight 
loops). In particular, this implies the same statistics for 
the L-PA and F-PA BCs on a long cylinder. In a broad 
strip (I > 7), the differences between L-PA and F-PA 
BCs must be related to the approach of the constrained 
path to the principal minimizer over a sequence of scales. 
Unlike local minimizers, principal minimizers are inde- 
pendent of the direction in which they are traversed. We 
expect that bulk segments of the L-PA curves are well 
described locally by SLE K . The conditioning of paths 
used in defining the DMP may be related to the proper- 
ties of the minimizing paths [23| . We also note that the 
LERW with reflecting BCs passes the same set of tests 
of conformal invariance as the L-PA 2DISG and fails the 
same set of tests of SLE K . 

In conclusion, we have numerically sampled over geo- 
metric objects in a system with disorder, domain walls in 



the 2D Ising spin glass, and tested their statistical geo- 
metric properties. We find that the domain walls pass to 
the left of a given point with probabilty consistent with 
SLE, wind around long cylinders in a manner consistent 
with conformal invariance, and that the sequences of con- 
formal maps that generate DWs, i.e., Loewner evolutions 
£(t), give a diffusion constant K — 2.32 ± 0.08 in accord 
with a fractal dimension df = 1.28 ± 0.01. We directly 
study the domain Markov property: it fails in small sys- 
tems, but we can not rule it out in larger systems. This 
set of tests, whose utility is validated by application to 
curves in LERW and MST, provides strong numerical 
support for a description of spin-glass domain walls with 
unconstrained endpoints by SLE, implying both confor- 
mal invariance and a domain Markov property on long 
scales. Domain walls starting from a localized bond are 
not consistent with the simplest form of SLE, though 
more complex conformally invariant descriptions, e.g., 
SLE with drift such as SLE K;P should be investigated. 
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